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CHAPTER  I 


INTRODUCTION 

An  integral  part  of  the  work  performed  on  tMs  contract  has 
been  the  development  of  a  spatially  resolving  electron  density  microwave 
diagnostic  technique  useful  in  a  range  of  electron  densities  infrequently 
explored  by  other  investigators.  This  technique  permits  accurate 
spatial  resolution  in  plasma  gradients  for  electron  densities ,  orders  of 
magnitude  below  that  possible,  using  a  conventional  microwave  inter¬ 
ferometer.  After  an  initial  calibration  phase  in  which  argon  was  the 
driven  gas  in  the  shock  tube,  this  microwave  technique  h.*s  been  used  to 
measure  electron  densities  and  electron  density  profiles  behind  normal 
air  shocks  in  the  range  of  Mach  number  8.5-11  and  initial  pressure 
conditions  01  ~  1  mm  Hg.  When  the  results  of  the  measurements  were 
compared  with  the  values  of  electron  densities  predicted  on  the  basis  of 
the  accepted  semi-analytical  rate  chemistry  program,  an  order  of 
magnitude  discrepancy  was  found.  Specifically,  the  measurements 
indicate  an  electron  density  behind  the  shock  which  tends  to  the 
equilibrium  value  at  a  rate  which  is  about  an  order  of  magnitude  faster 
than  the  one  predicted  by  the  numerical  calculations.  Since  the  shock 
conditions  where  the  measured  discrepancy  exists  are  relevant  to  the 
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150,000  ft.  altitude  range  aerodynamic  calculations,  it  is  deemed 
pertinent  to  explore  the  reasons  for  the  observed  discrepancy.  To  this 
effect  an  extensive  program  of  measurements  has  been  initiated  with 
different  gas  compositions  and  the  sources  of  errors  in  the  measuring 
technique  and  the  interpretation  of  the  results  have  been  analyzed.  The 
information  presented  in  this  report  will  not  necessarily  represent  the 
chronological  development  of  the  expe  imentation,  but  rather  a  summary 
of  the  physical  problem  and  the  methods  of  solutions. 

The  second  part  of  this  report  will  deal  with  the  properties 
of  an  electromagnetic  wave  propagating  through  a  plasma  with  random 
jr  periodic  spatial  density  variations.  Of  particular  interest  in  this 
problem  are  the  conditions  when  the  plasma  is  "under  danse"  and  the 
period  of  the  spatial  variation  is  about  equal  to  the  half  wavelength  or  a 
multiple  wavelength  of  the  incident  radiation.  The  experimental  research 
for  this  part  of  the  program  has  been  severely  handicapped  by  the 
inability  to  gererate  experimentally  a  plasma  which  has  the  required 
fluctuation  properties  and  is  amenable  to  an  analytical  evaluation.  Some 
tentative  results  of  this  work  will  be  presented. 
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CHAPTER  II 
THE  FACILITY 

SECT.  A- 1  INTRODUCTION 

The  explicit  purpose  of  the  following  detailed  description 
of  the  GASL  combined  shock  tube  and  waveguide  facility  is  to  enable 
other  investigators  to  duplicate  any  or  all  parts  of  this  highly 
successful  experimental  arrangement.  Since  it  is  the  details  rather 
than  the  basic  principles  that  are  time  consuming  in  the  fabrication 
of  such  a  facility,  great  effort  will  be  made  to  present  all  salient 
features. 
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SECT,  A- 2  BASIC  SHOCK  TUBE 

In  this  experimental  facility,  the  driven  section  of  the 
shock  tube  is  also  used  as  a  circular  waveguide.  Because  of  this 
duplicity,  some  precautions  must  be  taken  as  to  the  choice  of 
material,  the  surface  finish,  and  the  method  o'  connecting  tubing 
sections. 

A  compromise  material  for  the  tube  is  seamless  stainless 
steel  pipe.  While  this  material  is  certainly  not  the  best  electrical 
conductor  as  would  be  desirable  from  the  waveguide  point  of  view, 
it  offers  the  ad',  antages  of  strengh,  machinability, corrosion 
resistance  and  availability.  The  stock  material  used  for  the  driven 
sections  of  the  shock  tube  were  commercially  available,extra  heavy  . 
three  inch  stainless  steel  pipe.  This  tube  has  a  nominal  inside 
diameter  of  2.  9"  and  outside  diameter  of  3.  5". 

Highly  desirable  from  the  waveguide  aspect  is  that  the 
tube's  inner  diametei  be  constant.  For  use  as  a  shock  tuoe,  the 
inner  surface  should  be  smooth  with  no  wall  discontinuities  at  the 
tube  section  junctions.  To  satisfy  these  conditions,  each  ten  fool 
length  of  pipe  was  reamed  and  honed  by  a  specialist  in  this  field. 

The  resultant  inside  diameter  for  all  four  driven  sections  was 


measured  to  be  2,  917"  +  ,  001  ", 


The  problem  of  coupling  the  tube  sections  is  summarized 
as  follows:  for  use  as  a  waveguide^ood  and  complete  electrical 
contact  must  be  maintained  at  each  junction,  the  sections  should 
join  with  a  minimum  angle  between  their  axes,  discontinuities  in 
the  wal,  should  be  minimized,  and  the  junction  must  bo  pressure 
and  vacuum  tight  and  mechanically  strong  Also  desirable  is  a 
connection  that  can  be  easily  removed  and  reassembled,  since  the 
insertion  of  various  test  sections  may  be  required  at  any  junction. 

In  addition,  each  section,  should  he  able  to  be  rotated  and  be  locked 
independently.  This  requirement  originates  from  the  fact  that  the 
microwave  radiation  is  polarized  and  both  parallel  and  perpendicular 
measurements  are  sometimes  desirable. 

To  satisfy  these  requirements,  the  coupling  assembly 
shown  in  figure  A-2.  1  was  devised.  In  this  scheme,  one  end  of  each 
tube  section  is  machined  with  a  right  hand  thread  and  cylindrical 
surface  concentric  with  the  inner  axis.  On  the  other  end  of  each 
section, a  left  hand  thread,  an  identical  cylindrical  surface,  and  a 
partical  "O"  ring  groove  are  machined.  This  "O"  ring  forms  the 
seal  between  sections,  and  is  held  in  place  by  the  inner  surface  of 
.ne  coupling  nut.  These  are  relatively  simple  machining  operations. 
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and  if  done  in  the  order  statedjC.oncentr  icity  and  diametric  tolerances 
are  maintained  without  the  usual  distortion  problems  associated 
with  weld  type  coupling  structures. 

Bronze  is  used  for  the  coupling  nut  to  prevent  galling  of 
the  surfaces,  and  is  available  as  standard  unfinished  bronze  bearing 
stock.  If  the  tolerances  of  the  cylindrical  surfaces  are  held  to  .  002", 
the  miximum  angular  deviation  between  sections,  excluding  bending, 
is  approximately  .  1  degree.  Great  care  and  patience  should  be 
exercised  in  the  assembly  of  tube  sections  in  order  to  prevent  bend¬ 
ing. 

This  coupling  arrangement  satisfies  all  of  the  conditions 
previously  mentioned.  However,  in  order  to  assemble  and  dis¬ 
assemble  the  tubes,  the  downstream  sections  must  be  able  to  slide 
on  their  mount!  lg  pads.  Since  the  reaction  forces  on  the  shock  tube 
exert  only  axial  stresses,  there  is  no  need  to  yigidly  restrain  the 
tube  along  its  length,  so  long  as  the  main  reactionary  forces  are 
absorbed  at  the  diaphragm  location.  Since  this  is  not  only  feasible 
but  desirable  from  the  point  of  view  of  dimple  electrical  isolation 
from  ground,  the  supports  for  ‘he  driven  section  need  carry  only  • 
the  compressional  weight  load  of  the  tube.  This  being  the  case,  no 
extensive  foundations  or  structure  are  required.  Since  perfect 
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alignment  of  these  supports  to  the  accuracy  dictated  by  the  coupling 
is  not  practically  attainable,  an  adjustable  set  of  "vee"  jaws  can  be 
mounted  on  simple  "I"  beam  structures  to  support  the  entire  driven 
tube.  Such  a  sek  of  jaws  are  shown  in  Fig.  A-2.  2. 

These  jaws  allow  independent  adjustment  of  height  and 
horizontal  position.  'he  vertical  adjustment  is  achieved  by  rotating 
the  central  screw  which  has  right  and  left  hand  threads  operating 
though  the  respective  threads  in  the  "V"  jaws.  The  vertical  load 
on  the  jaws  is  transmitted  directly  to  the  mounting  plate  on  which 
they  slide.  Horizontal  adjustment  is  made  by  moving  the  angle 
bracket,  containing  the  center  slot  guide  bearing,  in  the  range 
provided  by  the  slots  in  the  bracket.  The  jaws  are  finally  clamped 
to  the  angle  bracket  which  is  also  slotted  in  the  horizontal  direction. 
Satisfactory  suppo  t  is  accomplished  with  two  sets  of  jaws  for  each 
ten  foot  section  of  tube.  With  this  basic  mounting  and  coupling 
arrangement  the  entire  driven  tube  can  be  assembled  in  one  day  with 
the  help  of  a  precision  level,  an  inexpensive  transit,  suitable  "vee" 
blocks  which  will  fit  on  the  cylindrically  machined  part  of  the  tube 
section,  and  a  little  patience  and  luck. 
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SECT,  A- 3  FAST  OPENING  MECHANICAL  VALVE 

General  Characteristics 

This  valve  has  been  developed  and  constructed  to  re¬ 
place  the  diaphragm  in  a  small  shock  tube. 

The  conditions  of  operation  are: 

Driver  pressure  -  100  to  500  psi 
Driver  temperature  -  up  to  20C°C 
Opening  area  -  7  sq.  inches 

Seal  •  vacuum  tight 

Additional  conditions  which  should  be  satisfied  are  re¬ 
mote  operation  and  minimum  of  adjustment  or  replacement  of  parts. 
The  adopted  solution  has  only  one  moving  part,  a  piston,  that 
provides  two  types  of  seal:  an  "O"  ring  seal  for  "acuum  and  a  sleeve 
valve  for  fast  action. 

Sequence  of  Operations 

This  sequence  is  shown  in  figure  A- 3.  1. 

(  a)  The  valve  is  closed.  Pressure  in  chamber  A  i  i 
equal  to  pressure  in  chamber  B  (  driver)  .  The 
resulting  force  pushes  the  piston  towards  the  seal. 
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(  b)  Venting:  Pressure  in  chaml  er  A  goes  down  towards 
atmospheric;  the  driver  pressure  remains  constant. 
The  resulting  force  pushing  the  piston  drops  and 
then  reverses. 

(c)  Vacuum  seal  opening:  When  the  piston  starts  to 

move  towards  the  left  the  force  pushing  it  increases 
suddenly  because  the  pressure  starts  acting  over 
the  surface  between  vacuum  seal  and  sleeve  seal. 

(  d)  Fast  opening:  When  the  sleeve  goes  out  of  the  driver 
chamber  the  gas  starts  to  flow  from  driver  to  driven. 
The  piston  has  been  accelerated  along  a  distance 
equal  to  the  length  of  the  sleeve. 

(  e)  Stopping:  The  piston  compresses  the  gas  left  in 

chamber  A,  stopping  and  oscillating  around  a  posi¬ 
tion  which  is  a  function  of  the  final  pressure  of  the 
tube. 

(  f)  Closing:  After  firing,the  driver  and  driven  are 
evacuated:  the  gas  left  in  chamber  A  expands  and 
closes  the  valve, 

(  g)  Charging:  The  gas  is  fed  first  into  chamb'  r  A, 

closing  the  valve  if  step  (  f)  fails  to  do  so,  and  then 
the  gas  goes  into  the  driver  through  a  check  valve. 


12 


The  two  main  facts  to  be  considered  in  the  action  of  the 

valve  are: 

(  1)  The  valve  behaves  as  an  unstable  system  near  the 

firing  pressure.  The  force  pushing  the  piston 

jumps  from  zero  to  a  maximum  value  in  a  very 

.4 

short  time  (<  10  sec).  This  fast  action  cannot 
be  obtained  by  controlling  the  flow  of  gas  with 
conventional  valves. 

(  2)  The  fast  seal:  Only  after  the  piston  receives  ite 
maximum  velocity  is  the  communication  between 
driver  and  driven  established.  For  a  given 
acceleration  of  the  piston  the  opening  time  goes 
down  with  the  square  root  of  the  length  of  the  sleeve. 
The  fast  seal  is  not  tight  and  the  possible  leakage 
between  the  opening  of  both  seals  establishes  a  compromise  between 
length  of  sleeve,  acceleration  and  quality  of  seal. 

Computation  of  the  Main  Characteristics 

The  calculations  are  very  easy  if  we  disregard  secondary 


effects. 
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The  main  parameters  are: 

SA  area  of  piston 
Sg  area  of  the  vacuum  tight  seal 
Sg  transversal  section  of  driven 
PA  pressure  in  the  chamber  A 
P_  pressure  in  the  driver 
W  weight  of  the  piston 
e  length  of  the  sleeve 

The  pressure  in  the  driven  is  ,  in  normal  working 
conditions,  in  the  order  of  a  few  millimeters  of  mercury  and  we 


disregard  it  in  the  calculations.  When  charged  it  is  PA  =  Pg,  and 
the  force  closing  the  seal  is  F  =  PB  SB’  The  Chamber  A  is  vented 
towards  atmospheric;  Pg  remains  constant  and  PA  goes  down 


reaching  the  firing  pressure  Pf  corresponding  to  zero  force  acting 

Sr 

on  the  pisi  Pj.  =  Pg  (  1  — - ) ;  if  Pf  is  equal  to  atmospheric 

SA  1 

pressure,  the  corresponding  PR  is  the  minimum  at  which  the  valve 


can  operate  for  a  given  rations., 

SA 


When  the  piston  starts  to  move,  the  force  pushing  it 


back  jumps  to  =  Pg  (  -  Sg)  -  P^  S^;  because  of  the  volume 

of  the  chamber  A  we  can  assume  that  PA=  Pf  during  the  fast  seal 
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opening  of  the  valve,  so  =  Pg  (  Sg  -  S^)  .  With  this  force 
we  can  compute  the  acceleration,  velocity  and  opening  time. 

We  consider  that  the  valve  is  open  when  the  passage 
area  of  the  gas  through  the  fast  seal  equals  the  transversal  area 
of  the  driven. 

The  corresponding  values  for  the  valve  in  operation  are: 
P  =  60  psi. 

Pg  =  300  psi. 

=  16  in2 
Sg  =  13  in2 


Fj  =  1  800  pounds 
opening  speed  =  120'/sec 
opening  time  =  300  microseconds 
Figure  A-3.  2  presents  the  calculated  shock  Mach 
numbers  as  a  function  of  the  pressure  ratio  between  driver  and 
driven  for  air.  This  calculation  assumes  an  infinitely  fast  opening 
diaphragm  and  no  energy  loss  mechanisms  occurring  in  the  flow. 
Experimental  points  (  +  )  have  been  superimposed  on  the  plot 
indicating  the  behavior  obtained  using  the  mechanical  valve.  The 


agreement  equals  or  exceeds  that  obtainable  using  a  bur  sting  diaphragm. 


CVCH3  X  1J3  D. VISIONS  PER  INCH 
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Modifications 

One  modification  refers  to  the  elimination  of  the 

minimum  firing  pressure.  The  section  S  is  equal  to  S0,  and  th_' 

A  B 

valve  remains  closed  after  the  total  venting  of  the  chamber  A. 

To  fire,  the  chamber  B  is  connected  to  the  area  Sg  -  Sg,  initiating 
the  motion;  in  this  way,  a  better  synchronization  is  obtained  and 
is  lower,  increasing  the  accelerating  force. 

The  gas  left  in  chamber  A  is  not  enough  to  stop  the 
piston,  so  additional  gas  is  injected  after  the  valve  is  open.  This 
is  achieved  with  a  bypass  conduit  which  the  piston  uncovers  on  its 


way  back  during  opening. 
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SECT.  A-4  ^HOCK  TUBE  PLUMBING 

The  governing  factor  in  the  engineering  of  the  shock  tube 
plumbing  is  saiety.  Since  the  driver  gas  normally  used  in  the  facility  is 
hydrogen,  the  possibility  of  explosion  either  internally  or  externally  due 
to  a  leak  is  omnipresent.  While  the  basic  tube  is  rated  for  over  1000  psi, 
the  various  detectors  and  fittings  are  in  general  not  designed  for  these 
pressures.  Hence,  even  a  limited  internal  explosion  is  serious. 

The  replacing  of  the  conventional  diaphragm  system  with  the 
mechanical  valve  has  considerably  reduced  the  hazards  of  operation.  With 
the  valve  the  use  oi'  the  shock  tube  has  become  a  closed  operation.  Since 
tube  need  not  be  opened,  the  only  air  entering  the  system  is  introduced 
during  the  -lontrolled  filling  cycle  o;  the  oriven.  This  air  charging  can  only 
be  performed  if  the  driver  is  not  charged,  the  driven  previously  evacuated, 
and  then  the  air  charge  cannot  exceed  on  atmosphere.  The  driver  cannot 
be  filled  until  the  pressure  in  the  driven  is  again  reduced  below  50  mm  Hg. 
Once  the  driver  has  been  filled  the  only  operafon  that  can  be  performed  is 
either  "fire"  or  "purge".  After  firing,  a  pressure  switch  in  the  driven 
starts  the  purging  cycle  automatically,  venting  excess  pressure  over 
atmospheric  to  the  atmosphere.  Then  the  vacuum  system  is  opened  and 
the  tube  evacuated.  While  it  is  possible  to  bypass  these  automatic  opera¬ 
tions  to  obtain  unusual  conditions,  it  is  generally  during  the  routine 
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operation  of  such  a  system  that  a  careless  error  is  serious.  Specific 
details  of  the  plumbing  system  are  not  warranted  here  due  to  the  fact  that 
the  use  of  the  mechanical  valve  instead  of  a  diaphragm  has  permitted 
additional  simplifica*  ions.  The  basic  system  to  be  installed  is  outlined 
below. 

The  combined  driver-mechanicai  valve  as  previously 
discussed  need  have  only  one  solenoid  type  filling  valve,  and  one  small 
mechanical  vacuum  valve  to  be  used  only  if  the  assembly  is  to  be  opened 
for  cleaning  or  servlvh  lg.  Another  small  solenoid  valve  may  be  used  if 
the  possibility  of  a  high  pressure  dump  is  desired.  The  driven  tube  is 
equipped  with  two  semi-automatic  valves.  The  first  'Is  an  over  pressure 
valve  whic''  will  seal  vacuum  but  wiU  open  a  two  inch  port  if  the  static 
pressure  in  the  driven  exceeds  one  atmosphere  absolute.  This  valve 
mechanically  interlocks  to  the  vacuum  valve  so  that  if  the  over  pressure 
valve  is  open,  the  vacuum  valve  must  be  closed.  The  vacuum  valve  is 
pneumatically  operated. 

With  thi  system  a  typical  operating  cycle  would  be  as  follows. 
The  driven  tv  be  is  pumped  down  by  opening  the  pneumatic  vacuum  valve.  The 
tube  is  filled  with  the  desired  test  gas  which  .s  limited  in  pressure  by  the 
over  pressure  valve.  The  tube  is  now  evacuated  to  the  desired  pressure. 

The  driver  is  then  filled,  It  may  be  noted  that  if  the  driver  had 
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been  inadvertently  tilled  and  fired  during  any  of  tlie  above  operations  all 
systems  including  the  vacuum  pump  would  have  been  protected  by  the  over 
pressure  valve.  When  the  firing  sequence  is  s  arted,  a  relay  is  locked 
so  that  the  next  pressure  rise  occurring  in  the  driven,  will  open  the  driven 
vacuum  valve.  After  firing,  if  the  pressure  exceeds  one  atmosphere,  the 
operation  of  the  vacuum  valve  is  blocked  until  a  safe  pressure  in  the 
driven  is  reached,  after  which  the  system  will  be  evacuated  leaving  it  in 
a  safe  condition  and  ready  for  the  next  test.  For  the  above  operation,  it 
was  assumed  that  the  main  mechanical  valve  operating  in  place  of  the 
diaphragm  closed  after  firing.  In  the  event  that  a  malfunction  occurs,  no 
further  operation  of  the  shock  tube  is  possiDle  until  the  trouble  is  correct¬ 
ed,  but  neither  are  any  safety  problems  present. 

Two  basic  static  pressure  gauges  have  been  found  to  be 
adequate  in  the  system.  To  monitor  the  driver  pressure  a  standard  250 
psi  Heiss  gauge  is  used  without  valves  or  restricting  orifices.  An  effective 
and  accurate  gauge  for  measurement  of  pre-firing  pressure  in  the  driven 
is  the  Alphatror  gauge.  This  gauge  (an  old  model)  has  three  linear 
ranges,  0  100  microns,  0-1  mm.  Hg.  ,  and  0-10  mm.  Hg.  ,  with  a  zero 
shift  raising  the  upper  range  to  20  mm.  Hg.  Supplementing  the  Alphatron 
gauge  is  a  standard  thermocouple  gauge  with  one  probe  in  the  driven  tube 
and  another  monitoring  the  press  are  at  the  vacuum  pump.  Neither  gauge 
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system  is  electrically  damaged  by  operating  at  atmospheric  pressure.  So 
far,  it  has  not  been  necessary  to  protect  any  gauge  systems  by  isolation 
during  firing.  This  is  important  from  a  safety  standpoint,  as  well  as  an 
aid  in  avoiding  possible  readout  errors. 
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SECT.  A- 5  VACUUM  SYSTEM 

A  satisfactory  range  of  initial  pressures  in  the  driven  tube  was 
deemed  to  be  from  50  microns  to  10  mm.  Hg.  To  avoid  the  additional 
complication  of  a  diffusion  pump,  a  high  quality,  high  capacity  mechanical 
pump  is  used.  Though  pumping  speed  is  slow  near  ultimate  pressure, 
this  range  is  seldom  used  for  decontaminating  the  tube  during  firing.  With 
the  tube  clean  enough  to  avoid  excessive  gassing  in  the  operating  pressure 
ranges,  the  simplest  ami  quickest  way  to  achieve  a  contamination-free  test 
gas  is  by  repeated  filling  with  test  gas  and  evacuating.  For  instance,  with 
the  tube  at  atmospheric  and  full  of  hydrogen  after  firing,  only  two  or  three 
minutes  pumping  will  evacuate  it  to  .  1  mm.  Hg.  By  refilling  with  test  gas 
and  evacuating  to  test  pressure,  remaining  hydrogen  contamination  is 
approximately  only  one  part  in  ten  thousand. 

With  the  addition  of  a  liquid  nitrogen  trap,  the  gassing 
contaminants  important  at  low  pressures  (  including  pump  oil)  are  re¬ 
moved.  However,  if  the  trap  is  used,  a  bypass  system  should  also  be 
employed  since  when  pumping  at  high  pressures,  severe  and  wasteful 
boiling  occurs  in  the  trap.  With  the  system  tight,  pumping  for  twelve  hours 
with  only  the  mechanical  pump  will  reduce  the  pressure  to  a  few  microns. 
Using  the  trap,  one  micron  can  be  obtained  in  two  hours. 
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JECT.  A-6  INSTRUMENTATION 

The  instrumentation  discussed  below  are  those  devices  and 
recording  instruments  pertinent  t  j  the  general  operation  of  the  shock  tube. 
The  microwave  instrumention  is  discussed  in  Section  A-7.  The  specific 
problem  under  study  requires  a  high  degree  of  accuracy  in  all  measuring 
and  recording  instruments.  Three  specific  areas  will  be  discussed; 
accurate  location  of  the  aerodynamic  shock,  measurement  of  dynamic 
pressures  and  recording  of  fast  rise,  small  amplitude  electrical  signals. 

Location  of  Shock  Front 

Accurate  estimates  of  the  equilibrium  properties  of  the  gas 
behind  the  incident  shock  can  be  made  if  the  initial  pressure  and  the  shock 
velocity  are  known,  the  latter  parameter  being  the  most  sensitive. 
Measurement  of  the  shock  velocity  can  be  made  by  measuring  the  time 
interval  between  two  points  along  the  tube  where  an  observable  property 
associated  with  the  shock  front  can  be  detected.  The  difficulty  in  this 
problem  is  the  measurement  of  a  property  directly  attributed  to  the  shock 
front.  Basically  only  three  fundamental  parameters  are  implicitly  as¬ 
sociated  with  the  shock;  a  temperature  change,  a  pressure  change  or  a 
density  change.  Measurement  of  parameters  induced,  such  as  electron 
density  and  light  radiation  to  locate  the  shock  front,  assumes  that  these 
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properties  abruptly  change  across  the  front.  This  assumption  is  never 
strictly  valid  but  may  be  satisfactory  at  high  Mach  numbers,  when  the  use 
of  ionization  gauges  or  photomultipliers  may  be  acceptable. 

If  it  is  desirable  to  detect  the  actual  aerodynamic  shock  front, 
one  of  the  fundamental  property  changes  must  be  observed.  A  direct 
measurement  of  the  density  change  can  be  made  optically  by  measuring 
the  deflection  of  a  light  beam.  However,  even  for  moderate  strength 
shocks,  this  method  is  very  delicate. 

A  direct  measurement  of  pressure  is  possible  except  that  the 
response  time  of  the  detector  strongly  limits  the  possible  resolutions. 

That  is,  using  a  good  quality  crystal  detector,  a  rise  time  of  three 
microseconds  is  typical.  If  the  shock  velocity  is  on  the  order  of  3  mm/ 
Usee  (~M9),  the  spatial  resolution  is  about  1  cm.  Assuming  the  rise 
time  is  consistent  for  two  gauges  and  the  shock  pressure  profile  is  constant 
during  the  interval  of  the  measurement,  a  velocity  measurement  is 
possible  with  a  spatial  resolution  of  .  5  cm.  Hence,  to  measure  velocity 
to  1%,  the  distance  between  gauges  should  be  .  5  meters.  It  is  clear 
however  that  for  the  purpose  of  locating  the  shock  as  a  method  of  cor¬ 
relating  events,  the  gauges  have  limitations.  In  addition,  if  measurements 
are  to  be  made  over  a  wide  range  of  conditions,  continual  adjustment  of 
sensitivity  must  be  made.  A  more  serious  d’fficulty  is  the  isolation  of 
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the  gauge  from  the  mechanical  shock  generated  when  the  diaphragm  breaks 
(valve  opens)  .  If  the  mechanical  shock  reaches  the  gauge  prior  to  the 
arrival  to  the  aerodynamic  shock,  a  spurious  signal  may  render  the 
measurement  useless.  A  discussion  of  pressure  gauges  and  measure¬ 
ments  will  be  made  in  a  following  subsection. 

Heat  Transfer  Gauges 

A  satisfactory  solution  to  a  detector  for  locating  the  shock 
front  has  been  found  to  be  a  heat  transfer  gauge.  The  basic  gauge  is  a 
thin  film  resistance  thermometer  element  bonded  to  a  glass  slug.  The 
gauge  response  is  a  function  of  not  only  the  temperature,  but  also  the 
density.  This  combined  response  allows  the  useful  operation  from  M~4 
to  ~14  in  this  facility.  The  gauge  is  fabricated  by  painting  a  thin  strip  of 
metallic  silver  coating  less  than  1  micron  thick  and  1  mm.  wide  on  a  glass 
disk  1/4  inch  in  diameter.  This  film  is  then  baked,  and  connecting  leads 
are  fed  through  two  holes  in  the  glass  disc  and  soldered  to  the  film.  The 
entire  front  surface  is  then  given  a  thin  coating  of  magnesium  flouride, 
which  aids  ir.  preventing  the  collection  of  the  charge  diffused  in  front  of 
the  shock  (Ref.  1)  .  Unsatisfactory  measurements  are  obtained  without 
the  coating.  In  operation,  a  steady  current  of  10-50  ma.  is  passed  through 
the  resistance  element  (thin  film).  The  temperature  change  as  the  shock 
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passes , creates  a  change  in  resistance  which  gererates  a  change  in  voltage 
if  the  current  is  nearly  constant. 

An  extremely  simple  method  of  generating  an  essentially 
constant  c’-rrent  is  to  use  a  large  voltage  with  a  large  series  resistor 


to  the  gauge  as  shown  in 
Fig,  A-6.  1.  To  be  sure, 
this  system  is  wasteful 
of  power,  but  its  simplicity 
over  a  suitable  transistor 
stabilized  constant  current 
supply  for  each  gauge  is 
obvious. 


Fig.  A-6. 1 
HTG  Power  Supply 


To  obtain  a  significant  pulse  (  20  volts)  from  this  circuit, 
very  large  amplification  (  20  x  10J)  of  the  voltage  pulse  must  follow.  If 
the  fast  response  of  the  gauge  is  to  be  fully  utilized,  the  amplifier  should 
have  a  rise  time  of  .  3  fi  sec . 

To  make  the  gauge  useful  as  a  standard  and  reliable  shock 
front  detector  the  circuit  shown  in  Fig.  A-6,  2  has  been  designed  and  built. 
The  circuit  is  intended  to  generate  a  sharp  pulse  ani  no  attempt  is  made 
to  preserve  the  long  term  (  10  y  sec.  or  more)  variations  occurring  at 


27 


the  gauge.  That  is,  all  the  time  constants  are  made  conveniently  small. 
The  grounded  grid  first  stage  utilizes  the  low  impedance  character  of  the 
gauge  to  simultaneously  provide  a  very  stable  first  amplifier  and  the 
required  bias  current  for  the  gauge  element.  The  following  two  stages 
of  amplification  are  conventional  RC  pentode  amplifiers  utilizing  an 
extremely  high  transconducUnce  tube.  The  output  cuthode  follower 
allows  extended  cabling  without  loss  of  response  time  which  is  typically 
.3  microsec.  The  combination  of  series  condensers  and  the  normally 
closed  relay  in  the  final  output  provides  a  simple  method  of  avoiding 
spurious  trigger  signals  which  occur  from  the  various  switching  opera¬ 
tions  performed  in  the  firing  sequence.  Because  of  poor  grounding 
practices  in  other  nearby  facilities,  a  common  failing  in  many 
laboratories,  minor  power  line  noises  generated  by  switching  transients 
are  amplified  by  the  high  gain  amplifiers.  To  minimize  this  problem, 
the  first  step  in  the  firing  sequence  disconnects  all  unnecessary  AC 
components  from  the  system  {fans,  heaters,  pumps,  etc).  The  first 
step  also  applies  a  DC  voltage  to  the  time  delay  circuit  of  the  output 
relay  which  opens  a  few  seconds  later.  This  system  is  not  a  replace¬ 
ment  for  proper  grounding  procedures  but  rather  an  insurance  to  avoid 
60  cycle  pickup  and  line  transients. 
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Care  must  be  taken  in  the  general  circuit  construction 
because  of  the  high  gain.  The  general  circuit  construction  utilizes 
Vector  turret  type  t  ockets.  Separate  filament  supplies  are  used  as  a 
precaution  against  high  frequency  oscillation.  The  separate  trans¬ 
formers  are  practical  if  multiple  units  are  to  be  constructed  on  the 
same  chassis.  For  one  amplifier  unit  one  supply  is  probably  adequate. 

The  B+  supply  can  be  unregulated  with  minor  filtering.  At  300  volts, 
the  current  consumption  is  about  70  ma.  per  amplifier.  No  gain  con¬ 
trol  is  used  since  in  general  the  purpose  is  to  provide  a  fast,  large 
trigger.  For  most  of  the  shock  conditions  shown  in  fig.  B-4.  5,  the 
amplifier  is  saturated,  thus  generating  a  60  voJt.  .  2  microsecond 
rise  time  pulse,  independent  of  shock  conditions. 

With  this  heat  transfer  gauge  and  amplifier,  spatial 
resolution  of  the  location  of  the  front  is  comparable  to  the  actual 
discontinuity  in  physical  parameters  occurring  across  the  shock. 

Typically,  this  resolution  is  1  mm,  indicating  that  for  a  shock  velocity 
of  3  mm/microsec.  the  required  distance  between  detectors  in  order  to 
determine  velocity  to  1%  is  about  10  cm,  assuming  the  timer  can  resolve 


.3  microsec  or  better. 
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Pressure  Detectors 

Measurement  of  pressures  or  pressure  variations  behind  the 
shock  is  extremely  difficult  if  accuracy  and  space  resolution  are  required. 
The  only  suitable  gau ge  elements  for  this  purpose  are  piezoelectric 
crystals.  The  basic  problem  in  the  design  of  a  gauge  for  shock  tut  t  use 
is  that  the  requirements  of  high  natural  resonant  frequencies  and  high 
sensitivities  are  contradictory.  The  second  problem  is  the  design  of  the 
crystal  mounting  in  the  gauge  so  that  the  output  is  insensitive  to  external 
stresses  at  right  angles  to  the  measured  stress. 

If  one  considers  that  the  useful  pressure  information  in  the 
shock  tube  occurs  for  an  interval  of  a  few  hundred  microseconds,  there 
'E  a  partial  elimination  of  the  measuring  problems  if  a  gauge  intended 
only  for  dynamic  purposes  is  used. 

A  precision  calibrated  crystal  microphone  with  a  natural 
resonant  frequency  of  about  200  KC  is  the  Massa  model  M-213.  This 
microphone  has  good  sensitivity  and  directional  resolution,  and  can 
/ithstauJ  severe  overloads.  A  typical  pressure  trace  using  this  micro¬ 
phone  -s  shown  in  figure  A-6.  3.  A  low  pass  filter  of  RC  constant  20 
micr oseconds  has  been  used  to  remove  the  initial  ringing  of  the  crystal, 
as  shown  in  fig.  A- 6.  4. 
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Upper  Trace:  Massa  Microphone 
with  no  filtering  -  gain  15  db 
5  volt/cm. 

Lower  Trace:  Electric  Field 
Detector  -  50  mv/cm. 

Sweep  Time:  20  (is/cm. 

Initial  Pressure:  3  mm  Hg  -  Air. 
Shock  Velocity:  3000  meters/sec. 

Figure  A-6.3 


Upper  Trace:  Heat  transfer  gauge 
current  10  ma,  5  mv/cm. 

Lower  Trace:  Massa  Microphone 
with  RC  =20  /is  filter  -  gain  20  db 
5  volts/cm. 

Sweep  Time:  50  /is/cm. 

Initial  Pressure:  2  mm  He  -  N 

O 


Figure  A-6.4 


Shock  Tube  Instrumentation,  Typical  Traces 
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Since  this  crystal  microphone  is  only  useful  for  dynamic 
pressure  measurements,  its  calibration  is  difficult.  In  general  the 
calibration  supplied  with  the  microphone  is  adequate.  If  doubtful,  the 
microphone  can  be  mounted  in  an  inflated  balloon,  the  pressure 
statically  measured  with  respect  to  atmosphere,,  and  the  balloon 
punctured.  The  resultant  step  pressure  change  supplies  a  positive 
calibration  point.  Since  the  gauge  and  amplifier  are  relatively 
inexpensive,  their  use  is  strongly  recommended  if  similar  dynamic 
measurements  are  required. 
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SECT.  A-7  MICROWAVE  INSTRUMENTATION 

Microwave  circuit 

The  circuit  is  as  shown  in  Figure  A-7.  1 ’includes  &  microwave 
signal  generator  (  2-4  kmc,  60  mw)  ,  a  coupling  device  delivering  this 
signal  to  the  shock  tube,  detectors  of  electric  and  magnetic  fields,  a 
recording  system,  consisting  of  oscilloscopes  and  photographic 
cameras,  marker  generators,  wave  shapers  and  delay  generators. 
Sometimes  a  counter  is  used  to  measure  the  time  between  markers 
and  consequently  the  speed  of  the  shock. 

The  microwave  signal  generator  is  a  reflex  klystron  with 
external  tuneable  cavity.  This  klystron  is  powered  by  a  power  supply 
which  has  the  possibility  to  change  reflector,  grid  and  beam  voltage. 

The  frequency  stability  of  this  generator  was  good  enough  for  our 
purpose,  but  the  measurements  being  made  were  particularly  sensitive 
to  changes  in  the  amplitude  of  the  signal.  Normally  the  klystron  is 
switched  on  and  off  and  the  maximum  tolerated  short  term  change  in  the 
power  output  should  be  below  2%.  Among  the  various  methods  tested, 
the  one  using  negative  pulses  in  the  grid  together  with  a  Zener-stabilized 
bias  power  supply  was  found  o  be  the  most  convenient. 
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Coupling  System 

The  output  of  the  klystron  is  fed  through  an  attenuator  to  a 
coaxial  cable,  an  isolator,  and  from  this  to  the  coaxial  waveguide 
transition.  The  shock  tube  equivalent  to  a  circulai  waveguide  operating 
in  the  TEy  mode.  The  coaxial  to  waveguide  transition  is  shown  in 
figure  A-7.  2.  Matching  is  made  by  changing  the  penetration  of  the 
antenna  and  the  position  of  the  short  circuit  piston.  The  additional 
short  circuit  stub  allows  for  adjustment  of  the  bandwidth  changing  its 
length.  The  shocktube  must  be  vacuum  tight,  and  this  is  achieved  by 
using  a  phenolic  tube  around  the  antenna  and  a  bellows  on  the  piston. 

Matching 

The  matching  of  the  coaxial  line  to  waveguide  is  very 
important.  All  the  measurements  are  based  on  the  fact  that  the 
amplitude  of  the  microwave  signal  coupled  into  the  shock  tube  is  constant. 
This  condition  is  difficult  to  obtain  because  the  plasma  produces  a  strong 
reflection;  standing  wave  ratio  of  twenty  is  not  unusual.  This  reflected 
signal  must  be  absorbed  without  reflection  in  the  isolator.  If  reflections 
take  place,  they  superimpose  on  the  incoming  signal  with  a  phase  that 
changes  as  the  plasms  moves.  This  •  ffect  is  equivalent  to  a  fluctuation 
of  the  input  signal.  The  standing  wave  pattern  then  becomes  a  function 
of  the  position  of  the  plasma  or  reflecting  surface  which  is  undesired. 
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Electromagnetic  detectors  mounted  on  the  waveguide  surface 
pick  up  the  sum  of  both  standing  wave  and  fluctuation  patterns,  both 
having  the  same  periodicity.  A  detector  moving  with  the  reflecting 
surface  picks  up  the  fluctuation  only.  The  phase  between  the  two 
patterns  may  take  any  value:  if  it  is  coincident  with  one,  for  instance 
E,  the  output  from  this  detector  will  give  a  higher  standing  wave  ratio 
than  the  output  from  the  other.  If  the  phase  angle  between  the  two 
patterns  has  any  value  the  resultant  pattern  is  not  symmetric.  In  the 
computation  of  the  plasma  properties,  it  is  as  Burned  that  the 
electromagnetic  energy  is  perfectly  coupled  to  the  waveguide,  hence 
perfect  symmetry  must  be  obtained  to  utilize  the  measurement. 

Several  methods  have  been  used  to  get  the  desired  matching 
but  the  only  one  ac.'  .  _..a  simple  enough  is  the  following.  A  short 
circuit  piston  is  displaced  along  the  shock  tube,  simulating  a  strong 
reflecting  plasma.  Moving  with  the  piston  is  a  detector  picking  up  the 
magnetic  field  near  the  wall,  where  it  is  maximum.  The  output  from 
this  detector  gives  the  fluctuation  produced  by  the  mismatch  only. 
Successive  adjustments  of  the  antenna  to  waveguide  transition  are  made 
attempting  to  eliminate  this  fluctuation, thus  giving  a  constant  output 
when  t\e  piston  moves. 
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In  order  to  reduce  the  frequency  sensitivity  of  the  antenna, 
the  signal  from  the  generator  is  altered  slightly.  The  antenna  short 
circuit  stub  is  then  adjusted  until  a  minimum  change  occurs  in  the  out¬ 
put  of  the  piston  detector  during  this  frequency  change. 

The  limit  of  this  system  of  matching  comes  mainly  from 
incomplete  contact  between  piston  and  wall.  A  typical  value  of  the 
fluctuation  after  the  matching  has  been  performed  is  2%, but  this 
corresponds  to  a  smaller  value  when  the  plasma  is  reflecting  the  signal 
because  of  the  significantly  better  reflecting  properties  of  the  metallic 
piston. 

Detectors 

The  electromagnetic  wave  propagated  inside  the  tube  has  a 
wavelength  Xg  which  is  a  function  of  the  frequency,  mode  of  propagation 
and  diameter  of  the  tube.  Some  of  the  radiating  energy  is  reflected 
from  the  plasma  and  the  ratio  of  reflected  to  incident  power  is  a  func¬ 
tion  of  electron  density,  collision  frequency  and  profile  of  the  plasma. 
The  combination  of  incident  and  reflected  waves  combines  to  form  a 

system  of  standing  waves  along  the  tube.  This  pattern  has  a  period  of 

X 

8  .  When  the  plasma  moves  the  pattern  also  moves  with  the  same 
2 

speed,  provided  the  properties  of  the  plasma  remain  constant. 
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A  device  detecting  the  electromagnetic  field  and  located  in 
a  fixed  position  detects  a  signal  whose  amplitude  goes  from  a  maximum 
value  equal  to  the  sum  of  incident  and  reflected  fields  to  a  minimum 
equal  to  the  difference  between  these  two.  For  a  plasma  velocity  v 
this  fluctuation  has  a  frequency  fm  = 

B 

Two  types  of  detectors  have  been  used,  one  for  the  electric 
and  the  other  for  the  magnetic  field  amplitude,  as  shown  in  figure 
A-7. 3.  The  electric  field  detector  is  located  in  the  region  where  it 
takes  its  maximum  value,  picking  up  the  component  parallel  to  the 
antenna  and  using  a  small,  flush  mounted  surface  in  order  to  avoid 
introducing  anything  inside  the  shocktube  which  would  disturb  the  flow. 
This  surface  is  connected  to  a  coaxial  resonant  circuit  that  matches 
impedances  to  the  crystal  detector  and  short  circuits  the  currents 
produced  by  the  flow  of  charges  from  the  plasma  to  the  tube  walls.  The 
magnetic  field  detector  ic  a  slot  excited  by  the  current  flowing  in  the 
wall.  Because  of  the  relation  between  magnetic  field  and  current  this 
slot  is  equivalent  to  a  loop  located  in  a  plane  perpendicular  to  the  slot. 
The  field  to  be  checked  is  the  transversal  component  of  H.  The  slot 
should  be  oriented  in  a  transverse  section  of  the  tube  and  at  the  point 
where  the  electric  field  is  maximum.  The  signal  induced  in  the  slot 
is  coupled  to  a  two  line  resonant  circuit  and  to  a  crystal  detector. 
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For  both  E  and  H  crystal  detectors  an  output  is  obtained 

which  is  relate  o  the  amplitude  of  the  respective  field.  The 

characteristics  of  the  crystals  cannot  be  considered  linear  nor  quadratic 

and  the  required  accuracy  demands  a  calibration.  The  time  response 

-7 

of  the  detectors  is  less  than  10  sec  if  their  output  is  directly  connected 
to  cathode  followers.  For  simplicity's  sake,  AC  amplifiers  are  used 
exclucively  hence  the  DC  level  is  not  present.  However,  switching 
off  the  klystron  gives  this  zero  level  reference, permitting  the  absolute 
amplitude  determination  from  the  single  trace.  The  markers  are 
equally  spaced  at  known  intervals  providing  a  timing  reference  and  a 
means  of  compensating  for  nonlinearitieB  in  the  Bweep. 

To  perform  the  calculation  of  the  placma  properties  the 
electric  and  magnetic  fields  should  be  determined  in  the  same 
transversal  section  of  the  tube.  A  special  section  ha-  been  constructed 
allowing  the  placing  of  the  two  detectors  plus  a  shock  front  detector 
in  the  same  transversal  section.  In  this  arrangement,  when  the  signal 
from  one  detector  is  maximum  the  other  is  minimum  and  vice  versa. 

The  oscilloscopes  are  triggered  using  various  shock  front  detectors 
as  discussed  in  section  A-6. 
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SECT.  A-8  DATA  RECORDING 

The  only  economically  feasible  method  of  recording 
analogue  data  consistent  with  shock  tube  time  is  the  oscilloscope, 
although  it  is  possible  to  record  some  data  with  video  type  tape 
recorders.  However,  the  cost  and  resolution  favor  the  oscilloscope 
and  camera.  While  this  basic  technique  is  common-place,  to  achieve 
optimum  fidelity  of  the  data  consistent  with  the  errors  in  the  system, 
each  step  of  the  process  must  be  carefully  analyzed.  This  errcr 
analysis  "s  in  many  instances  a  function  of  the  recording  time.  The 
range  of  interest  for  his  facUity  are  sweep  times  of  from  5  microsec. 
to  500  microsec. 

Two  types  of  distortion  can  be  delineated  in  the  over-all 
process;  these  are  amplitude  and  time  distortion.  To  trace  the  source 
of  these  and  estimate  their  relative  and  absolute  magnitudes,  the  over¬ 
all  recording  process  is  divided  into  the  following  categories:  the  input 
amplifier  up  co  the  vertical  deflection  plates  of  the  CRT,  the  sweep 
genera,  or  and  horizontal  amplifier  up  to  the  horizontal  deflection  plates, 
and  the  distortion  in  the  CRT.  Also,  the  photographic  transfer  of  the 
CRT  face,  and  the  subsequent  extraction  of  numerical  data  from  the 
photographic  traces  must  be  considered.  3y  far,  the  greatest  offender 
in  these  distortion  processes  is  that  due  co  the  CRT  itself.  While  the 
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other  factors  are  not  always  unimportant,  each  of  the  other  parts  of 
the  orocess  except  numerical  data  reduction  can  be  maintained  to  an 
order  of  magnitude  better  than  the  CRT  distortion,  if  care  is  taken. 

Even  with  the  best  selected  tube;*,  in  order  to  obtain  distortion  levels 
on  the  order  cf  .  5%  only  the  central  region  of  the  CRT  screen  can 
be  utilized.  Of  course  by  so  doing, this  also  tends  to  reduce  all  the 
other  distortion  causes  listed.  However,  if  these  processes  are 
independently  investigated  it  becomes  clear  that  the  C;  T  is  the  arch 
criminal. 

The  specific  procedure  used  to  obtain  the  data  for  the  shock 
electron  density  profile  analysis  discussed  in  section  B-3  is  outlined 
below.  The  critical  traces  were  made  using  selected  CRT's  so  as  to 
minimize  distortion.  The  photographic  transfer  of  the  trace  was  made 
using  a  modified  scope  camera  to  reduce  the  image  size,  Polaroid 
10,  000  speed  film,  extremely  fast,  enables  the  use  of  a  small  lens 
aperture.  With  this  arrangement  no  distortion  of  a  precision  grid 
placed  on  the  CRT  was  observable.  The  picture  obtained  was  converted 
to  a  transparency  by  using  a  very  small  aperture  on  the  enlarging 
len6.  The  subsequent  transparency  was  enlarged  in  a  3"  x  4"  slide 
projector.  In  this  projector,  a  3/8"  aperture  was  placed  inside  the 
standard  lens  assembly  and  a  green  filter  was  added  to  reduce  chromatic 
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abberation.  The  resultan  olargement  of  about  12  diameters  was 
negligibly  distorted.  The  projection  screen  consisted  of  a  layer  of  fine 
frosted  glass  backed  by  Plexiglass.  The  Plexiglass  was  precision 
scribed  with  .  1"  squares.  The  scribed  lines  were  painted  with 
phosphorescent  paint.  By  using  a  "black  lite  "  fluorescent  behind  the 
screen,  the  scribed  marks  were  clearly  visible,  and  direct  transfer 
of  numerical  data  from  the  picture  could  be  taken  by  working  from 
behind  the  screen.  Such  traces  have  also  been  made  on  large  graph 
paper,  but  such  paper  is  very  often  distorted,  introducing  significant 
errors.  Under  the  best  conditions,  the  overall  error  from  amplifier 
input  to  data  output  was  estimated  to  be  .  4%  for  the  maximum  value 5 
increasing  to  2%  for  values  1/20  of  maximum.  The  increase  in  error 
becomes  significantly  large  for  small  numerical  values  because  the 


trace  cannot  be  resolved. 


CHAPTER  III  ELECTRON  DENSITY  MEASUREMENT 
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SECT.  B  -  1  ANALYTICAL  FORMULATIONS  FOR  DETERMINING 
ELECTRON  DENSITY 

The  following  is  a  summary  of  the  analytical  expressions 
used  to  determine  the  electron  density  in  the  combined  shock  tube- 
waveguide.  The  polarized  radiation  is  introduced  ;n  the  downstream 
er.d  of  the  tube.  Probes  to  detect  the  radial  electric  field  and  the 
tangential  magnetic  field  are  upstream  of  the  antenna,  as  described 
in  Section  A-7.  To  determine  the  point  by  point  variations  in  electron 
density  as  the  plasma  moves  past  the  detector  station,  expressions  are 
required  relating  the  measured  field  quantities  to  the  local  electron 
density. 


For  the  TE^  mode  of  propagation,  the  radial  electric 
field  and  the  tangential  magnetic  field  at  the  surface  of  the  waveguide 
are  given  by, 

iU)t 

Er  =  A,  J(y  ,R)  cos  0u(x)  e 


Hfl  *  — — 
6  imjo 


a'(x)  F 

u(x)  r 


(1) 

(2) 
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where  y  for  the  TE^  mode  is  — —  ,  R  being  the  tube  radius. 


The  function  u(x)  is  an  undefined  function  which ,  however, 


contains  the  plasma  properties.  It  must  be  a  solution  to, 
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u"  +  ulri_  k  -y  =  0 


where  primes  indicate  space  derivatives,  and  CO  is  the  radiating  frequency, 


with 


K  =  1  + 


TS) 


COp  and  V  being  the  plasma  frequency  and  collision  frequency  respectively. 
Assuming  the  field  components  are  measured  at  right  angles,  with  the  Er 
field  in  the  direction  of  polarization  of  the  radiation,  then 


Er  ■  AQu(x)  e 


_  u'(x) 


HU  - - rt  E  * -  u'(x)  e 

o  fj,  co  u(x)  r  jito 


Since  the  values  of  the  fields  that  are  measured,  assuming  a  square  law 
detection,  are  the  square  of  the  real  parts  of  the  respective  field,  the 


measured  quantities  are, 


2  sjc 

ErEr  ■  A0  u(x)  u  (x) 


*W  —Tz*™ 

M  co 


where  *  indicates  the  complex  conjugate.  Solving  for  in  Eq,  (3), 
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(9) 


The  solution  o i  in  terms  of  the  field  values  of  Eqs.(7)  and 

(8)  can  be  shown  to  be, 

:rEr  1 
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and 


Im  K2  = 


2  c^  2 
Re  K  -  yL 
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The  plasma  properties  in  terms  of  Re  K  and  Im  K  are  given  as. 


WP  2 

—  =  (1-Re  K2) 

or 


f  w  |  j 


and , 


(12) 


V_  _  Im  K2 
40  Re  KZ-1 


(13) 


The  values  of  the  field  in  the  above  equations  are  in  MKS  units. 
However,  it  may  be  noted  that  each  field  quantity  appears  in  a  ratio.  All 
that  is  required  for  solution  in  terms  of  the  nondimensional  measured 
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quantities  are  i.ie  appropriate  scaling  ratios. 

2  2 

If  the  measured  values  of  the  field  are  designated  as,  E  and  H  , 

(E^)*  (E2)" 

then  the  distance  axis  scaling  needed  for  the  ratios  - and  - = —  can 

E^  E^ 

be  determined  from  the  ratio  of  the  measured  wavelength  Xg  (in  units  of 
the  elongated  trace).  That  is,  if  the  horizontal  (distance  axis)  coordinate 
is  measured  in  units  of  Z  and  the  real  coordinate  in  X,  units,  then, 


hi  ma 


where  X  is  the  actual  waveguide  wavelength  (MKS  units).  Hence  to  scale 

g 

the  absolute  ratios  of  field  values  to  the  measured  quantities, 
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To  scale  the  combined  ratios  of  electric  and  magnetic  fields, 
it  may  be  noted  that  in  the  free  space  region,  the  ratios  of  the  maximum 
field  values  (occurring  180°  apart)  is  given  by, 


(H8H8*>max  __4lL 

,EW>ma*  " 
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If  the  ratio  of  the  same  measured  quantities  is  designated  by  (a),  i.e.  , 

,2 


H 


max 

,2  "  a* 

max 


then, 
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Then  Eqs.  (10)  and  (11)  can  be  put  in  terms  of  measured  quantities  to  give. 
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Hence,  the  plasma  properties  can  be  obtained  through  the  use  of  Eqs.  (12) 
and  (13), 

In  summary,  to  calculate  the  plasma  properties  at  any  point 
behind  the  shock  the  following  parameters  are  required. 


waveguide  wavelength 


H 


a  * 


max 


max 


S 


Xmhrn 

A  R 
co*  2rf 


taken  in  the  free  apace  region 


distance  scaling  ratios 

measured  wavelength  on  enlarged  trace) 

O 

R  is  tube  radius 
f  radiating  frequency 


and  E^,  (E^)',  (E^)",  the  field  values  and  derivatives  at  the 


wall  of  the  waveguide. 
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SECT.  B-2  DATA  PROCESSING  PROGRAM 

In  the  preceding  section  describing  the  analytical  formulation 
necessary  for  the  electron  density  measurement,  the  basic  quantities 
required  for  solution  are  the  real  value  of  electric  field  squared  (EE  ), 
its  first  and  second  derivatives,  the  real  value  of  magnetic  field  squared 
(HH ' ) ,  and  its  first  derivative.  It  was  shown  that  no  absolute  measure¬ 


ment  of  these  quantities  was  required,  if  the  ratio 


EE 


max 


could  be 


HH 


max 


obtained  in  the  free  space  region  in  front  of  the  shock.  The  precision  of 
the  final  calculation  is  extremely  sensitive  to  these  measured  quantities. 
While  it  has  been  indicated  in  Section  A-8  that  the  relative  accuracy  for 
part  of  the  recorded  data  is  about  .4%,  it  is  important  to  note  that  this 
accuracy  is  related  to  the  data  taken  at  that  instant  of  time.  Because  all 
of  the  components  of  the  measuring  system  are  constantly  varying ,  it  is 
essential  that  all  of  the  critical  information  be  obtained  from  one  set  of 
traces  taken  at  one  time. 

The  greatest  variables  in  this  system  are  the  crystal  detectors, 

although  sweep  rates  and  amplifier  gains  may  also  change.  If  the  desired 

% 

signal  to  be  detected  is  for  instance  EE  ,  the  crystal  detector  must  have 
a  square  law  behavior.  Since  this  requirement  is  never  precisely  true  a 


correction  must  be  made  on  the  output  signal.  It  will  be  shown  how  this 
correction  can  be  made  directly  from  the  numerical  data  derived  from  the 
single  trace,  thus  avoiding  the  need  for  corroborating  data  taken  at  a 
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different  time. 

The  calculation  of  electron  density  is  divided  into  three  oper¬ 
ations.  First,  the  raw  numerical  data  taken  from  the  scope  trace  is 
corrected  for  the  non-square  law  behavior  of  the  crystal.  Second,  from 
this  corrected  data,  first  and  second  derivatives  are  computed.  Third, 
the  ;cted  data  and  derivatives  are  combined  point  by  point  as  pre- 

scriuc^  oy  the  analytic  equations  of  the  preceding  section. 

Processing  of  the  Raw  Data 

The  basic  input  for  the  calculation  is  approximately  150  data 
points  from  each  of  the  raw  electric  and  magnetic  field  deLector  output 
traces.  These  traces  include  one  full  period  of  the  sinusoidal -like 
function  preceding  the  plasma. 

Because  careless  errors  maybe  introduced  in  the  transcription 
or  punching  of  the  raw  data,  the  program  first  checks  if  successive  points 
are  monotonically  increasing  or  decreasing.  The  object  of  this  procedure 
is  to  identify  data  that  is  clearly  impossible,  and  allow  corrections  or 
changes. 

Once  the  data  passes  this  initial  test  the  points  are  smoothed  by 
a  least  square  fitting  procedure.  In  this  computation,  eight  points  are 
used  to  obtain  the  constants  of  a  second  power  equation  by  least  squares  fit. 
Having  computed  the  equation  constants,  the  values  of  the  function  for  the 


center  two  points  are  calculated  and  stored  as  smooth  data.  The  program 
then  jumps  two  points  and  recomputes  the  least  square  tit  for  the  next 
eight  points  which  of  course  includes  six  from  the  previous  calculation. 
After  this  procedure  has  been  performed  for  all  data  points,  the  initial 
raw  data  is  destroyed  and  the  smoothed  values  inserted  in  its  respective 
place.  The  need  for  smoothing  the  data  is  comparable  to  the  use  of  a 
"French"  curve  in  graph  plotting.  The  basic  assumption  is  that  the  data 
is  originally  smooth  and  that  conversion  to  digits  has  introduced  a  random 
error  in  the  least  significant  place.  By  reducing  these  errors  an  improve¬ 
ment  can  be  obtained  in  the  subsequent  derivatives  that  must  be  taken. 

The  next  step  is  to  modify  the  data  for  the  incorrect  response 
of  the  detector.  This  correction  procedure  is  based  on  the  fact  that  the 
data  function  outside  the  plasma  must  be  a  sinusoid  if  the  detectors  are 
measuring  EE*.  The  basic  problem  is  to  find  a  multiplicative  function 
Cg(x)]  such  that  when  multiplied  by  the  data  function  [f(x)]  the  resulting 
function  will  be  a  sinusoid  elevated  above  zero. 

i.e.  ,  g(x)  f(x)  =  A  sin  0  +  A  +  Ot 

The  difficulty  in  determining  the  g(x)  function  is  that  since  a,  which  is 
related  to  the  standing  wave  ratio,  is  unknown,  then  there  are  an  infinite 
set  of  g(x)  functions  which  will  satisfy  the  above  criteria.  The  solution  to 
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this  dilemma  is  found  in  the  fact  that  the  crystal  nonlinearity  (to  the  square 
law)  must  be  a  smooth  and  slowly  varying  function,  a  fact  clear  from  the 
physics  of  the  semiconductor.  If  Ot  is  on  the  order  of  A/3  or  smaller,  the 
latter  condition  becomes  a  very  stringent  limit  on  possible  g(x)  functions. 
The  fi"st  step  to  determine  the  g(x)  function  is  to  determine  the  period  of 
the  sinusoidal  like  function  f(x)  in  the  free  space  region.  This  is 
accomplished  by  computing  a  function  like  a  derivative  over  the  entire 
curve.  This  derivative  procedure  consists  of  computing  multiple  differ¬ 
ences  at  each  point  and  averaging  the  result.  The  r-.ason  for  this  technique 
is  to  establish  the  zero  slope  points  with  the  least  error  being  introduced 
from  the  random  data  scatter. 

After  locating  the  first  two  zeros,  forty  points  are  generated 
by  interpolation  between  zeros,  which  of  course  correspond  to  the  first 
maximum  and  minimum  or  the  reverse  of  the  original  curve.  The  forty 
points  are  chosen  such  that  if  the  function  were  a  true  sine,  then  the 
increment  of  the  function  between  points  would  be  equal.  This  is  done  so 
that  the  g(x)  function  will  not  be  unduly  weighted  by  a  preponderance  of 
points  at  each  end.  The  first  trial  g(x)  function  is  generated  by  assuming 
the  perfect  sine  curve  has  the  same  magnitude  as  the  data  function.  Hence 
the  g(x)  function  at  maximum  and  minimum  (A  and  B  in  Fig.  B-2.1)  is 
identically  one.  This  then  produces  a  g(x)  typically  like  that  in  Fig.  B-2.2. 
On  the  g(x)  curve,  the  region  from  A  to  D,  where  D  is  arbitrarily  1/2  f(x) , 
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second  order  polynomial.  Point  C,  th“  value  of 
compute!  from  the  curve  fit  constants  and  this  v^lue 


is  least  square  fitted  by  a 
g(xl  for  f{x)  min.  is 
is  used  to  compute  the  value  of  ft  . 

Distorted  Sine  Function 


g(x)  Determination 


Using  thit  value,  a  new  magnitude  A  is  computed  and  the  entire 
g(x)  function  for  each  point  of  f{x)  is  generated.  This  second  trial  g(x) 
function  is  again  curve  fit  to  a  second  order  polynomial  by  method  of  least 
squares.  From  the  constants  of  the  polynomial,  the  value  of  g(o)  and  the 
standing  wave  ratio  are  calculated  and  the  difference  between  the  curve  fit 
g(x)  and  the  trial  computed  g(x)  points  is  printed.  At  this  time  in  the  pro¬ 
ceeding,  the  computer  is  tired  and  returns  control  to  the  operator.  The 
operator,  making  a  human  type  judgment  on  the  basis  of  the  differences  in 
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the  last  points  may  choose  a  different  value  by  a  few  per  cent  of  g(o)  and 
the  above  process  will  be  repeated. 

It  is  desirable  that  the  companion  curve,  HH  ,  to  the  EE  curve 
be  processed  up  to  the  determination  of  the  g(x)  correction  function  as 
indicated  above.  If  the  g(x)  functions  have  been  properly  chosen  for  each 
curve,  the  standing  wave  ratio  will  be  identical.  With  a  little  practice  the 
standing  wave  ratios  can  be  made  to  agree  to  better  than  2%. 

Once  the  g(x)  functions  are  satisfactorily  determined,  all  the 

$  $ 

data  is  corrected  thus  yielding  che  "true"  values  of  EE  and  HK  along  the 
curves.  By  a  method  of  multiple  differences,  the  derivative  of  these 
functions  are  taken  yielding  (EE*)',  (EE*)"  and  (HH*)'. 

The  computer  program  thus  far  describee  is  intended  to  take 
data  derived  from  a  single  trace,  apply  corrections  and  produce  derivatives 
at  fixed  intervals.  The  functions  generated  must  not  only  be  self-consistent, 
but  also  should  be  consistent  with  the  data  derived  from  the  companion 
curve.  (EE"  and  HH  form  a  set}  It  has  been  shown  how  the  program 
"forces"  an  overall  amplitude  correction  by  assuming  the  detected  function 
in  the  free  space  r  ;gion  is  sinusoidal.  This  type  of  correction  will  com¬ 
pensate  for  most  amplitude  distortions  due  to  any  cause  that  is  not  a  function 
of  the  time  coordinate  (horizontal  position  on  the  CRT). 

Additional  information  which  can  be  obtained  from  the  data  is 
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the  relative  wavelength  of  the  sinusoidal  function  in  the  free  space  region. 

It  may  be  replied  that  the  wavelength  determination  was  necessary  to 
establish  the  period  of  the  "forcing"  sinus  function.  Since  the  waveguide 
wavelength  of  the  microwave  radiation  can  be  accurately  measured  prior 
to  testing,  it  is  used  as  the  cnly  additional  piece  of  information  that  is  not 
directly  obtained  from  the  picture.  If  the  measured  wavelength  and  the 
relative  wavelength  are  combined,  ?i  overall  scaling  is  achieved.  In 
addition,  if  the  scaling  factors  for  the  EE*  and  HH*  are  compared,  an 
accurate  contraction  or  expansion  factor  can  be  obtained  and  used  to  provide 
exact  time  correlation  between  points  on  either  curve.  If  one  additional 
marker  (such  as  the  location  of  th*  shock  front)  is  common  to  both  curv  s 
an  exact  relation  between  the  curves  can  be  made.  If  the  trace  also  contains 
time  markers,  the  shock  velocity  can  be  obtained.  It  may  be  noted  that  the 
time  markers  also  provide  the  location  of  the  zero  amplitude  axis  at  the 
instant  the  trace  was  made,  a  procedure  which  removes  the  uncertainty  due 
to  drift  in  all  systems. 

An  additional  piece  of  information  required  for  the  calculation 
of  electron  density  is  the  ratio  of  the  maximum  on  the  electric  field  trace 
in  the  free  space  region  to  the  maximum  on  the  magnetic  detector  trace 

EE*rnax 

Hir*  .  This  ratio  provides  all  the  scaling  information  needed  to  use  the 

max 

nondimensional  curve  data  in  an  absolute  determination  of  electron  density. 
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It  may  be  recalled  that  the  maximum  values  required  were  obtained  to  pro¬ 
vide  amplitude  information  for  the  "forcing"  sinusoidal  correction  function. 

Hence,  it  has  been  demonstrated  that  the  complete  set  of 
information  required  to  compute  the  electron  density,  with  ;he  exception  of 
the  waveguide  wavelength  is  directly  obtainable  from  the  two  detector 
curves.  The  advantage  of  this  procedure  is  that  no  absolute  measurements 
are  required,  with  the  exception  of  the  waveguide  wavelength  and  interval 
between  time  markers.  (The  latter  is  used  only  for  the  velocity  deter¬ 
mination.)  This  clearly  implies  that  no  absolute  calibration  and  maintenance 
of  that  calibration  is  needed  for  the  instruments  used  in  recording.  Except 
for  the  requirement  of  a  linear  time  scale  (assuming  shock  velocity  constant), 
an  overall  correction  of  most  of  the  distortions  occurring  in  all  phases  of 
the  data  collection  process  is  obtainable  from  the  data. 

All  of  the  information  obtained  in  the  "data  processing"  phase 
of  the  program  is  subsequently  used  in  the  point  by  point  determination  of 
electron  density  as  per  the  formulations  of  the  preceding  section.  Since 
this  computation  is  entirely  straightforward  no  discussion  will  be  given 


here. 
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SECT.  B -3  EXPERIMENTAL  DETERMINATION  OF  THE  ELECTRON 
DENSITY  PCHIND  A  NORMAL  AIR  SHOCK 

The  following  is  a  brief  description  of  the  experimental  conditions 
for  which  the  longitudinal  electron  density  profile  his  been  determined.  A 
series  of  six  independent  sets  of  data  have  been  processed  corresponding  to 
nearly  identical  shock  conditions.  The  purpose  of  the  nearly  redundant  data 
was  to  check  the  overall  repeatability  and  the  effect  of  small  variations  in 
shock  condition  about  the  mean. 

The  specific  shock  conditions  are  tabulated  below. 


Case 

Shock  V  elocity 
(meter/sec) 

Initial  Pressure 
(mm  Hg) 

Note 

1 

3100 

2.5 

Upper  curve 

2 

2900 

1.5 

Lower  curve 

3 

3030 

2.0 

4 

3015 

2.05 

Mean  values 

5 

3035 

1.95 

plotted 

6 

3020 

2.0 

The  velocities  quoted  above  were  calculated  directly  from  the  electric  and 
magnetic  field  detector  traces  as  previously  described  in  section  B-2.  The 
microwave  instrumentation  used  is  described  in  Section  A-7.  The  specific 
waveguide  wavelength  for  all  of  the  runs  was  20.2  cm. 

One  typical  set  of  traces  is  shown  in  Fig.  B-3,1.  Note  the 
periodic  time  markers  which  also  give  the  aer  o  level  of  the  signals.  The 
location  of  the  shock  front  was  determined  from  a  third  trace. 


Set  of  Electric  and  Magnetic  Field  traces  used  as  -starting 
data  in  the  calculation  of  the  electron  density  profile  behind  the 
shock.  (Case  4  in  table). 


Figure  B -3. 1 
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The  numerical  data  from  these  curves  was  obtained  in  the  manner 
described  in  Section  A-8.  Subsequent  numerical  analysis  was  carried  out 
using  che  computer  program  described  in  Section  B-2.  The  results  for  the 
six  cases  are  summarized  in  Fig.  B-3,2.  Here  the  mean  curve  represents 
the  results  for  cases  3  through  6  which  correspond  to  the  same  velocity 
conditions  ±  .  3%.  Cases  1  and  2  correspond  to  higher  and  lower  velocities 
respectively.  It  may  be  noted  that  the  behavior  of  the  electron  density  is  in 
correspondence  with  these  conditions.  The  values  of  equilibrium  electron 
density  were  calculated  for  these  shock  conditions  using  the  procedure 
outlined  in  B-4  and  Fig.  B-4.2. 

Because  the  shock  front  marker  detector  used  in  these  runs  was 
derived  from  a  pressure  detector  as  described  in  Section  A-6,  r  small 
d'  placement  uncertainty  existed  in  the  results  consistent  with  the  "jitter" 
of  the  relatively  slow  rise  of  the  detector.  The  "mean"  curve  in  Fig.  B-3.2 
represents  a  repositioning  of  the  separate  results,  where  the  maximum 
shift  corresponded  to  ,75  cm. 

The  test  gas  for  all  of  these  cases  was  commercial  compressed 
breathing  air.  The  basic  result  has  not  been  significantly  modified  by  any 
of  the  following  procedures. 

(1)  Addition  of  a  liquid  nitrogen  "freczeout"  trap  on  the  test 
gas.  Such  a  trap  would  of  course  remove  all  water  and  oil  present. 
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(2)  Liquid  nitrogen  trap  in  the  vacuum  system  to  remove  any 
trace  of  pump  oil  from  the  tube. 

(3)  Mixing  of  high  purity  oxygen  and  nitrogen  in  proportion  to 
"normal"  air  and  slight  variation  about  the  standard  mixture.  These  gases 
were  also  liquid  nitrogen  trapped  to  remove  water  or  oil. 

(4)  The  design  of  the  mechanical  valve  discussed  in  Section  A-3 
was  motivated  by  the  desire  to  remove  the  possibility  of  contamination  due 
to  the  use  of  plastic  diaphragms.  However,  no  significant  change  of  the 
electron  density  profile  was  observed  after  the  conversion  to  the  mechanical 
value . 

Thus  far,  there  is  no  indication  of  any  type  that  the  result 
obtained  is  not  indicative  of  the  electron  density  profile  for  the  given  shock 
conditions.  Ultimately,  the  results  from  the  spectroscopic  analysis  of 
the  plasma  region  behind  the  shock  as  discussed  in  Section  B-7  will  deter¬ 
mine  if  a  significant  impurity  exists. 


SECT.  B-4  SUPPORTING  INFORMATION 


Radial  Uniformity  of  the  Plasma 

One  assumption  tr»->de  in  the  analysis  of  Section  B-l  for  the 
electron  density  was  that  of  radial  uniformity  of  the  plasma  density. 

Such  uniformity  can  be  altered  by  two  factors.  First  is  the  presence  of  a 
boundary  layer  along  the  walls  in  the  region  behind  the  <?hock.  However, 
for  the  penetration  distances  indicated  by  the  results,  this  should  have  an 
insignificant  effect  on  the  assumption  of  uniformity. 

The  second  effect  is  that  of  a  ^ .  ipolar  diffusion.  A  calculation 
of  this  effect  for  two  operating  conditions  has  been  performed  as  sum¬ 
marized  in  the  following  subsection.  The  result  indicates  the  assu  nption 
of  radial  uniformity  is  basically  correct. 


Analytical  Determination  of  Transverse  Electron  Density  Profiles 

The  ambipolar  diffusion  equation  for  electrons  in  cylindrical 
coordinates  is: 
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where  D  =  — -  is  the  ambipolar  diffusion.  No  $  deper.  .ence  of  the 

electron  distribution  is  included. 


Only  diffusion  normal  to  the  wall  will  be  studied  and  for  that 


particular  case  Eq.  (1)  becomes’ 
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I  JL 

r  dr 


(2) 


The  term,  ne  which  is  related  to  the  creation  of  electrons  due  to  chemical 
reactions  inside  the  ionized  gas,  as  yet  must  be  defined. 

The  mechanism  for  the  creation  of  the  electron  is  given  by  the 
motion  of  a  shockwave  along  a  tube.  From  the  several  reactions  which 
can  generate  electrons  in  air  due  to  the  effect  of  the  shockwave,  the  most 
important  is: 

N  +  O  -  NO+  +  e  (3) 


Therefore , 

n  =  k  (K  ru.n^-n  2) 
e  r'  c  N  O  e 


(4) 


Replacing  (4)  into  (2)  we  have: 
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Eq.  (5)  can  be  written  in  the  form: 
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where neQ  is  the  equilibrium  electror  number  density. 

It  is  known  from  chemistry  considerations  (see  Fenner, 
Chemistry  Problems  in  Jet  Propulsion),  that, 
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Then  Eq.  (6)  becomes: 


I  _1 

r  dr 


k  n 
r  eo 

4.  . 

D 


d-y2) 


o 


(8) 


This  equation  has  a  boundary  condition  given  by  the  following  relation  (see 
GASL  TR-316,  Chapter  II). 


ne  dr 


(9) 


which  can  be  written. 


D_  <±£ 
y  dr 


(10) 


The  Eq.  (8)  with  the  condition  (10)  must  be  solved  numerically  and  it  has 
been  done  for  5  shock  wave  configurations.  The  data,  corresponding  to 
these  configurations  appear  in  Table  I. 


TABLF  I 


Shock 


Velocity 

T 

dp 

Ko 

neo 

Of 

m/sec 

°K 

-3 

cm 

cm 

2750 

2790 

8.45  x  10“3 

4.4  x 

109 

1.404 

3000 

3000 

1.00  x  10“2 

2.10  x 

1010 

.639 

3250 

3210 

1.08  x  10“2 

7.2  x 

10 

.348 

.500 

3390 

1.22  x  10*2 

2  x 

lO11 

.204 

3750 

3660 

1.30  x  10‘2 

6.2  x 

10n 

.120 

Of  primary  interest  is  the  effect  of  the  coefficient. 


i  k  n 
1  r  eo 

&  "  D 


(1 


upon  the  electron  distribution. 

The  coefficient  Of  has  the  dimension  of  X  .length  and  can  be 
considered  to  be  a  characteristic  diffusion  length.  Table  I  shews  that  Of 
changes  greatly;  it  implying  that  for  the  cases  where  Of  is  small,  the 
diffusion  will  be  very  slow  at  the  beginning  but  will  fall  off  sharply  near 
the  wall  of  the  cylinder.  This  does  not  happen  for  much  higher  values  of  Of 
(see  Fig.B-4.1).  The  experimental  investigation  made  at  GA3L  for  a  case 
where  the  shock  wave  velocity  was  3100  m/sec  corresponds  to  large 
values  of  Of  .  For  this  velocity  case,  the  transverse  gradients  are  small 
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over  most  of  the  tube  radius,  becoming  large  only  in  a  thin  outer  sheath. 


Sensitive  Check  on  the  Equilibrium  Electron  Density 

By  assuming  a  simple  model  for  the  plasma  behind  the  shock  to 
be  a  sharp  transition  from  free  space  to  equilibrium  electron  density  levels, 
a  formulation  for  the  standing  wave  ratio  due  to  the  addition  of  the  trans¬ 
mitted  and  reflected  signals  from  such  a  plasma  can  be  easily  made. 

By  combining  the  formulation  of  the  S.W.R.  with  the  properties 
of  the  gas  behind  an  air  shock  as  a  function  of  shock  velocity  and  initial 
pressure,  a  comparison  with  experimental  results  can  be  made.  This 
procedure  is  summarized  below. 

First,  the  standing  wave  ratio  will  be  calculated  in  terms  of 
the  electron  density  and  collision  frequency  which  are  assumed  uniform  in 
the  plasma  region. 

The  generalized  propagation  constant  as  defined  by 


is  given  by, 


d2E  2 
— -  +  PTE  =  0 

dx^ 


;  “p 


K  = 


■> 

-  v 
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where,  40  it:  the  propagating  frequency 
40p  the  plasma  frequency 
U  the  collision  frequency,  and 

y  is  the  root  of  the  Bessel  function  for  the  TEj^  circulator 
waveguide  mode  divided  by  the  tube  radius. 

The  propagation  constant  can  be  formed  as  follows.  Define 


K  =  k  K 
°  g 


where 


k 


o 


Now  let  the  radial  electric  field  and  the  tangential  magnetic  field  in  the 
free  space  region  be  given  by 


-ik  x  -ik  x 
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and  in  the  uniform  plasma  region, 


E 
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Normalizing  A  -  1,  the  boundary  conditions  are  matched  yielding, 


-ik  x  ik  x 

E=e  +  0  e 

r 


with 


0  = 


1  -  K 
1+K 


SL 

g 


Thus,  the  real  component  of  the  field  produced  is, 


*  *  *  "ik«x  ik  x 

E  E  "  1  +  00  +  00  e  e 
r  r 


This  it  reduced  to, 

ErEr*  -  [l±/(Re0)Z+{Im0)2]2 

The  two  solutions  give  the  maximum  and  minimum  values.  Thus  the 
standing  wave  ratio  is , 

r  / - 12 

1+y,Re6)?  +dm0)2 

SWR  =  v  — 

1  -/V*2 +IIm6>2 

where  0  is  a  function  of  ajp  and  v  . 

By  solving  the  normal  set  of  aerodynamic  equations  and  using 


a  programmed  Mollier  diagram,  the  conditions  behind  the  shock  can  be 


given  as  a  function  of  shock  velocity  and  initial  pressure. 

To  simplify  the  calculation  of  equilibrium  electron  density 
from  temperature  and  density,  standard  tables  were  programmed  f  the 
computer.  Thus  as  a  function  of  shock  velocity,  the  equilibrium  electron 
density  is  computed  with  the  initial  pressure  as  a  parameter.  Figure 
B-4.2  shows  this  plot.  The  collision  frequency  is  computed  from  the 
approximation 

5x  1012p 

V  =  - 

J~T 

with  p  the  pressure  in  atmospheres,  T  in  °K. 

If  now  these  plasma  properties  ar  e  used  to  compute  K  ,  the 

8 

standing  wave  ratio  can  be  calculated  using  the  previous  formulation  as 
a  function  of  sho  ;k  velocity  with  the  initial  pressure  as  a  parameter 
(Fig.  B-4.3).  Hence  as  a  function  of  accurately  determinable  quantities, 
the  electromagnetic  characteristics  of  the  entire  system  can  be  computed 
and  subsequently  compared  with  the  experimental  quantities. 

It  should  be  appreciated  that  such  a  comparison  is  a  very 
sensitive  check  on  the  assumptions  made.  First,  the  aerodynamic 
properties  behind  the  shock  must  precisely  follow  the  gas  dynamic  theory. 
Second,  the  electron  density  and  collision  frequency  must,  on  the  average, 
be  identical  with  the  equilibrium  values.  Third,  the  assumed  geometric 
uniformity  of  the  plasma  region  must  be  nearly  correct.  Fourth,  the 


4-4- H 


calculation  of  the  SWR  based  on  the  simple  step  model  should  be  right. 
When  comparing  with  experimental s ,  the  variables  measured  must  be 
identically  the  same  as  those  used  in  the  formulation.  For  instance, 
the  electric  field  detector  after  correction  should  indeed  be  measuring 
only  EE*. 

In  this  light,  the  correspondence  between  the  analytically 
determined  and  the  measured  SWR  as  presented  in  Fig.  B-4, 4  is 
remarkable.  Note  here  that  since  the  driver  pressure  is  constant  for  all 
these  tests,  in  order  to  change  the  shock  velocity  there  is  a  corresponding 
cnange  in  initial  pressure.  Hence,  Fig.  B-4  4  corresponds  to  a  cross 
plot  of  Fig.  B-4.  3  using  the  operating  characteristic  curve  in  Fig.  B-4.  5. 
While  by  itself  this  agreement  provides  probable  verification  of  the 
assumptions  made,  in  conjunction  with  the  electron  density  profile  meas¬ 
urement  discussed  in  Section  B-3,  the  measurement  provides  strong  proof 
that  the  point  by  point  technique  is  indeed  measuring  the  ionization 
relaxation  length  behind  the  shock. 
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SECT.  B-5  RATE  CHEMISTRY  PROGRAM  TO  DETERMINE  PROPERTIES 
BEHIND  A  SHOCK 

The  following  is  a  brief  description  of  a  computer  program  intended 
to  determine  the  properties  of  the  flow  downstream  from  a  point  where  all 
states  are  specified.  Since  the  parameters  which  determine  the  downstream 
properties  are  functions  of  many  variables  the  complete  set  of  equations 
governing  this  system  must  be  selved.  If  the  parameters  are  "rate  dependent", 
then  the  set  must  be  sol/ed  by  a  step  by  step  procedure.  The  complexity  of 
this  system  may  be  appreciated  by  realizing  that  at  each  step,  the  chemical 
reactions  of  the  species  (up  to  50  reactions  and  20  species)  must  be 
numerically  accounted  for,  where  the  corcentrations  of  the  species  are  under 
the  simultaneous  effect  of  other  chemical  reactions,  departures  from  vibra¬ 
tional  equilibrium,  and  variations  in  prescribed  properties  of  the  flow. 

The  validity  of  the  result  obtained  downstream  depends  on 
correctness  of  the  assumptions  made  for  the  various  dynamic  processes 
involved  and  the  numerical  constants  used  to  describe  the  energy  states  of 
the  species  and  the  rates  of  reactions. 

For  the  most  part,  the  assumptions  used  have  been  satisfactorily 
verified  in  other  physical  regimes.  Likewise,  the  energy  states  for  each 
species  are  calculated  from  "fundamental"  data,  /here  the  quotation  mark 
means  order  of  magnitude  of  better  accuracy  than  the  error  in  the  calculations. 
The  third  set  of  parameters,  the  dynamic  chemical  rate  constants,  are. 


unfortunately,  poorly  documented. 

The  difficulty  in  determining  these  rate  constants  is  that  even  for 
a  simple  chemical  system,  a  number  of  reactions  can  be  occurring.  At 
the  temperatures  and  pressures  of  interest,  the  relaxation  times  are 
generally  in  the  low  microsecond  region.  Confounding  this  difficulty  is 
that  for  some  species,  for  example  NO,  the  reaction  constants  must  be 
established  in  a  complicated  environment,  where  some  knowledge  of  the 
processes  and  rates  of  the  other  species  must  be  assumed. 

This  program  has  been  carefully  checked  for  internal  consistency 
and  has  been  checked  for  the  shock  conditions  for  which  the  important  rate 
constants  were  deduced.  The  results  for  the  latter  demonstration  case 
are  shown  in  Fig.  B-5.1,  The  relaxation  distance  for  this  case  corre¬ 
sponds  to  a  fraction  of  a  centimeter  behind  the  shock. 

The  reactions  considered  and  the  rate  constants  used  are  given 


in  Appendix  I. 
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SECT.  B-6  COMPARISON  BETWEEN  EXPERIMENTAL  AND  COMPUTED 
ELECTRON  DENSITY  PROFILE  BEHIND  MACH  8.9  SHOCK 

The  basic  shock  conditions  for  which  this  comparison  has  been 
made  were  detailed  in  Section  B-3.  The  specific  velocity  and  initial 
pressure  correspond  to  Mach  8.9  shock  (3025  m/s)  at  initial  pressure  of 
134,000  teet  altitude  (2.0  mm  Hg) .  The  curve  in  Fig.  B-6.1  was  taken 
from  the  "mean"  determination  discussed  in  Section  B-3. 2. 

The  analytically  determined  curve  in  Fig.  B-6.1  was  calculated 
from  a  program  similar  to  that  discussed  in  Section  B-5.  Actually  three 
separate  programs  using  slightly  different  procedures  and  constants  have 
yielded  essentially  the  same  curve.  The  obvious  order  of  magnitude 
difference  in  Fig.  B-6.1  warrants  serious  consideration. 

Experimentally  the  checks  discussed  in  Section  B-3  have  been 
performed  to  test  the  validity  of  the  result.  Basically,  if  an  error  exists 
in  the  experimental  information  it  is  either  due  to  the  fact  that  the  electron 
density  measuring  technique  is  incorrect  or  that  the  plasma  is  not  air 
alone,  but  has  a  contaminant. 

The  additional  measurements  discussed  in  Section  B-4  demonstrate 
quite  conclusively  that  the  equilibrium  electron  density,  assuming  pure 
air  is  used,  is  achieved  in  a  short  distance  compared  to  the  waveguide  wave¬ 
length  of  the  probing  microwave  frequency.  This  result  is  entirely  consistent 
with  t.ie  point  by  point  measurement.  In  addition,  since  '.he  calculation 
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of  Section  B-4  was  made  assuming  a  normal  air  mixture,  if  an  impurity 
were  present,  the  only  reasonable  effect  it  could  have  would  be  to  act  as 
a  catalyst.  The  rational  for  this  conclusion  is  that  since  tue  difference 
between  the  measured  and  calculated  SWR  in  Fig.  B-4.4  corresponds  to 
about  20%  in  electron  density,  which  ic  well  inside  the  expected  error,  any 
contaminant  present  does  not  add  significantly  to  the  equilibrium  density. 

If  this  is  so,  the  contaminant  must  either  act  as  a  catalyst  or,  by  a 
remarkable  coincidence,  have  the  same  equilibrium  density  as  air. 

If,  on  the  other  hand,  the  error  is  in  the  calculated  electron 
density  profile,  such  an  error  must  be  effective  for  the  temperature  and 
density  conditions  pertinent  to  this  shock  condition,  since  the  programmed 
calculation  gives  satisfactory  agreement  with  experiment  for  much  higher 
shock  Mach  numbers. 

The  obvious  location  of  such  an  error  is  in  the  rate  constants. 
Since  it  has  only  recently  become  possible  to  alter  the  rate  constants  and 
expression  used  in  the  program,  no  detailed  discussion  of  this  work  is  given 
here.  However,  it  maybe  noted  that  if  an  error  exists  in  low  energy 
reaction,  a  compensating  error  may  be  present  in  a  high  energy  reaction. 
Thus,  at  high  temperature  the  first  error  maybe  undiscernable  since  the 
result  is  principally  governed  by  the  higher  energy  reaction.  However,  at 
lower  temperatures  (recalling  that  the  rate  expressions  are  exponential  in 
nature),  the  error  in  the  low  energy  reaction  becomes  important. 
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Presented  in  Fig.  B-6.2  are  various  experimentally  determined 
values  for  the  forward  rate  constant  for  the  as3umed  principle  ionization 
process,  NO  -♦  NO+  +  e.  On  this  figure,  a  range  of  values  for  this  reaction 
calculated  from  the  properties  of  the  measured  profile  are  given.  The 
range  indicated  is  due  to  the  uncertainty  in  the  other  processes.  No  strong 
significance  should  be  placed  on.  this  initial  estimate  aside  from  the  fact 
that  it  is  not  inconsistent  with  the  other  measured  points. 

Two  basic  areas  of  work  are  suggested  by  this  result.  One  is 
the  re-evaluation  of  the  chemical  rate  processes  so  as  to  be  consistent 
with  the  high  and  low  Mach  number  conditions.  With  the  new  flexible  com¬ 
puter  program  empirical  estimated  changes  can  easily  be  made.  Also,  a 
check  for  an  impurity  in  the  air  plasma  should  be  made.  Such  a  check  can 
be  attempted  spectroscopically  as  discussed  in  the  next  section,  number  B-7. 
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SECT.  B-7  SPECTROSCOPIC  SEARCH  FOR  IMPURITY 

In  the  previous  section,  it  was  suggested  that  one  possible 
"error"  in  the  electron  density  profile  measurement  was  that  an  impurity 
existed  in  the  shock  gas  which  made  the  result  not  representative  of  normal 
air.  Although  all  possible  attempts  have  beer  made  to  minimize  impurities 
(see  Section  B-3)  a  persistent  criticism  of  *'ae  measurement  is  the  possible 
existence  of  an  impurity. 

To  investigate  this  possibility,  a  grating  spectrometer  of  large 
effective  aperture  has  been  obtained.  (Jerrel  Ash  model  75-000.)  At  the 
time  of  writing  of  this  report,  initial  runs  are  being  made  using  the  shock 
tube. 

Extensive  calibration  of  the  instrument  has  been  made  using  an 
rf -excited  discharge  tube.  This  discharge  tube  was  designed  so  that  any 
given  gas  or  impurity  could  be  introduced  and  tested  at  known  pressures. 
While  the  conditions  of  this  source  are  not  the  same  as  the  shock  tube,  the 
basic  energy  levels  of  the  species  are,  of  course,  unaltered,  even  if  the 
relative  intensities  are  dissimilar.  Hence,  if  a  specific  impurity  is 
suspected,  extensive  tests  can  be  made  to  establish  its  detailed  spectra 
which  can  then  be  used  to  verify  or  eliminate  the  existence  of  such  an 
impurity  in  the  shocked  gas  spectra. 

A  subordinate  piece  of  equipment  to  the  spectrograph  is  a 


recording  densitometer.  A  somewhat  unique  and  inexpensive  densitometer 
has  been  built  using  the  basic  parts  of  the  precision  enlarging  facility  used 
to  obtain  numerical  data  from  the  scope  trace. 

The  spectrographic  plate  (5x7)  is  loaded  in  a  modified  plate 
holder  for  the  3x4  projector.  One -half  of  the  seven-inch  plate  is  enlarged 
about  15  diameters.  This  enlargement  reduces  the  requirement  for  high 
precision  stages  and  slits  usually  required  in  such  an  instrument.  At  the 
normal  image  focal  plane ,  a  set  of  silicon  photo  cells,  one  with  a  slit,  is 
mounted  on  a  precision  platform  driven  by  a  synchronous  motor.  At  the 
lens  of  the  projector,  a  light  chopper  is  used  to  achieve  about  20J  cpk 
fluctuation.  The  use  of  the  fluctuating  light  source  allows  the  use  of  narrow 
band  ac  amplifiers  to  increase  the  output  from  the  photo  cells,  thus 
simplifying  the  design  of  the  amplifiers  and  minimizing  the  effects  of 
changes  in  ambient  lighting  conditions. 

Two  photo  cells  are  used,  one  with  a  narrow  slit  (.010")  to 
sample  the  trace,  and  the  other  with  no  slit  (l/4r  sq.  in  area)  to  sample 
the  unexposed  part  of  the  plate.  The  outputs  from  the  amplifiers  for  these 
two  photo  cells  are  put  in  bucking  series  in  a  diode  detecting  circuit.  The 
purpose  of  the  two  photo  cells  is  to  compensate  for  changes  i:i  "fogging" 
conditions  and  remove  zero  drift  due  to  changes  in  light  intensity.  In 
actual  use  no  discernable  drift  has  been  observed  even  for  intensity  changes 


of  a  factor  of  tw  >.  The  detected  compensated  output  is  recorded  on  a  10" 

synchronous  driven  strip  recorder. 

The  overall  system  from  spectrograph  to  output  chart  is 

repeatable  to  better  than  .1%  enabling  an  absolute  line  location  to  better 
o  0 

than  3.0  A  out  of  3000  A,  which  is  twice  the  minimum  resolution  of  the 


mtire  system. 
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CHAPTER  IV 

PERIODICALLY  DISTURBED  PLASMA 


SECT.  C-l  INTRODUCTION 


The  simple  criterion  used  to  determine  radar  cross 

section  of  an  ionized  wake  is  to  assume  that  fcr  electron  density 
(ji)  2  03  2 

regions  where  -E-<  1,  no  reflection  occurs,  and  if — ?—>  1  total 
03z  U3Z 

reflection  occurs.  This  simple  criteria  is  useful  if  the  collision 
frequency  is  low  and  tv  ~  density  gradients  are  large  through  the 
above  critical  conditions. 


The  problem  of  interest  here  is  the  electromagnetic 

03p2 

characteristics  of  a  plasma  region  where  — —  «  l,  but  where 

03 2  03p2 

stationary  or  time  varying  regions  of  higher  density  (——still  less 

032 


than  1)  occur  with  spacing  comparable  to  the  half  wavelength  or 
multiple  of  the  incident  radiation.  This  situation  is  comparable  with 
a  multilayered  dielectric  filter,  except  that  in  general  the  locations 
of  the  maximum  election  density  regipns  are  rardomly  changing. 

It  must  be  noted  that  while  it  is  possible  to  obtain  a  strong 
reflection  from  such  a  periodic  dielectric  stru:ture,  it  is  equally 
probable  that  total  absorption  will  occur,  by  the  same  process.  How¬ 
ever,  if  the  structure  is  randomly  varying,  the  process  of  interest  is 
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that  of  reflection,  which  ultimately  is  detected, 
energy  is  equivalent  to  a  transparent  region  as 
characteristics  are  concerned. 


The  totally  absorbed 
far  as  the  reflected 
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SECT.  C-2  EXPERIMENTAL  OBSTACLE 

In  order  to  properly  formulate  the  analytical  model  required 
to  generalize  the  problem,  me  basic  experimental  information  must 
be  supplied.  Of  course  such  information  must  be  sufficiently 
quantitative  that  a  true  model  can  be  formulated.  Therein  lies  the 
experimental  problem,  for,  since  the  basic  plasma  required  is  strongly 

f“p|* 

underdense.l — £-1  «  1,  measurement  of  the  plasma  properties 

I  u)  / 

becomes  extremely  difficult. 

One  approach  to  this  measurement  problem  is  to  induce 
the  required  perturbation  in  such  a  manner  thin,  the  plasma  properties 
can  be  accurately  deduced  from  other  parameters. 

One  such  experimental  arrangement  considered  was  to 
propagate  a  strong  sound  wave  into  the  shocked  gas.  With  such  an 
arrangement,  the  amplitude  and  spacing  of  the  periodic  structure  could 
easily  be  changed  within  limits.  Since  the  maximum  expected  varia¬ 
tions  in  pressure  would  be  small  compared  to  the  base  pressure,  the 
electron  density  variation  might  be  assumed  to  be  directly  proportional 
to  the  pressure  variation.  That  is 


e 
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While  all  parts  of  this  scheme  are  correct  in  concept, 
severe  problems  arise  in  coupling  the  required  sound  energy  to  the 
shock  tube.  Although  several  dilierent  coupling  systems  have  been 
tried,  the  attempt  to  use  sound  waves  to  induce  plasma  variations 
has  been  abandoned  because  of  the  extensive  facility  that  would  be 
required. 

To  obtain  an  order  of  magnitude  estimate  of  the  effect  of 
a  turbulent  wake  on  the  reflected  radiation,  an  experiment  was 
performed  in  the  shock  tube.  A  cylindrical  bar  was  mounted  across 
the  tube  parallel  to  the  direction  of  polarization  of  the  radiation. 
Several  diameters  were  tried  up  to  1/2",  which  was  the  maximum 
diameter  consistent  with  an  unchoked  flow.  Underdense  plasma 
cc editions  were  run  and  observations  were  made  of  the  reflected 
signal  with  the  use  of  che  electric  field  detector,  in  the  polarized  and 
antipolarized  directions. 

Even  though  the  scale  of  the  turbulence  expected  was  an 
order  of  magnitude  smaller  than  the  wavelength,  if  very  strong  effects 
were  to  be  associated  with  this  type  of  wake,  a  noise  signal  would  be 
imposed  on  the  reflected  power  signal.  Since  no  such  signals  were 
observed,  it  may  be  concluded  that  small  scale  effects  are  unimportant 


in  this  process. 


92 


SECT,  C-3  TENTATIVE  RESULTS 

An  attempt  was  made  to  see  if  strong  effects  could  be 
obtained  under  conditions  where  t  \e  scale  of  the  disturbance  was 
comparable  to  the  wavelength,  even  if  no  precise  description  of  this 
disturbance  could  be  given. 

The  basic  scheme  used  was  the  addition  of  a  "Tee"  section 
to  the  shock  tube,  where  the  area  of  the  arm  of  the  "Tee"  is  about 
1/2  the  area  of  the  tube.  As  the  flow  passes  the  "Tee"  section, 
under  certain  conditions,  a  "whistle "  type  effect  can  occur.  If  this 
is  the  case,  periodic  pressure  fluctuations  should  be  induced  in 
the  flow.  Such  pressure  variations  should,  in  turn,  create  regions 
of  higher  electron  density. 

In  figure  C-3.  1  a  trace  is  presented  showing  the  pressure 
fluctuations  occurring  in  the  arm  of  the  "Tee"  near  the  main  flow. 
Although  such  variations  could  only  be  obtained  for  a  very  restricted 
set  of  operating  conditions,  figure  C-3.  1  clearly  shows  that  some 
type  of  periodic  pressure  producing  mechanism  is  occurring. 

To  check  the  validity  of  the  result,  a  simple  step  model 
for  the  plasma  behind  the  shock  can  be  used  to  compute  the  reflected 
signal.  Assuming  the  pressure  defect  is  carried  with  the  flow,  the 
simple  model  for  the  plasma  postulates  a  reflection  from  the  front 
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of  the  shock  and  a  second  reflection  from  the  pressure  disturbed 
region  which  exists  behind  the  shock  a  distance  (z)  where  z=  vQt 
with  vQ  the  particle  velocity  with  respect  to  the  shock  and  time  t 
stai'ting  from  the  creation  of  the  disturbance. 

For  the  moment,  if  it  is  assumed  that  the  waveguide 

wavelength  is  not  appreciably  altered  in  the  first  plasma  region 

directly  behind  the  shock,  then  it  is  clear  that  when  the  distance 
X 

z  =  — -S-,  the  reflected  signal  from  the  disturbed  region  will  be  180 
out  of  phase  from  the  reflection  occurring  at  the  front  of  the  shock, 
thus  producing  a  minimum  in  the  standing  wa  »  ratio.  On  the  other 
hand,  for  z  =  Xg/2  the  two  reflected  signals  are  in  phase  and  a 
maximum  should  occur  in  the  SWR. 

In  figure  C-3.  2  a  trace  of  the  electric  field  is  shown  for 
the  critical  whistle  conditions  described.  While  this  result  must  be 
considered  tentative ,  it  is  clear  that  before  the  shock  enters  the  "Tee" 
section,  the  SWR  is  small  indicating  an  "underdense"  plasma.  After 
the  transition  of  the  "Tee"  where  the  effects  cannot  be  clearly 
described  at  the  present  '  ne,  the  amplitude  modulation  of  the  SWR 
is  seen  to  decrease  and  then  subsequently  to  increase.  The  end 
transition  occurs  as  first  the  piasma  passes  the  detector  and  then 
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Top  Trace:  Heat  transfer  gauge 
located  at  entrance  to  "Tee" 
section. 

Lower  Trace:  Pressure  in  arm 
of  "Tee". 

Sweep  Time:  100  fis/cm. 


Figure  C  -3.  1 


. .  «  *  f 
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Electric  Field  Detector.  Located 
3  meters  downstream  of  "Tee" 
section. 

Sweep  Time:  200  (Js/cm. 

Note:  The  cancellation  and 
subsequent  growth  of  the  signal  as 
the  "filter"  structure  goes  through 


condition. 


Figure  C-3,2 
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reaches  the  antenna,  where  the  signal  drops  to  zero.  Since  this  type 
of  modulating  behavior  is  consistent  with  the  simple  model  and 
inconsistent  with  other  possible  mechanisms,  ’t  is  felt  that  this  trace 
represents  the  t-'pe  of  strong  reflections  that  can  be  obtained. 
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APPENDIX  I 


SUMMARY  OF  CHEMICAL  REACTIONS  AND  THEIR  RATE  CONSTANTS1 


This  appendix  summarizes  a  large  number  of  chemical  reactions 

which  are  significant  at  the  high  temperatures  of  reentry  into  the  eaith's 

atmosphere  at  hypersonic  speeds.  In  all  cases  the  data  was  taken  from 

2 

the  survey  of  Steiger,  and  the  original  sources  may  be  found  in  this 
reference. 

The  presentation  here  includes  two  classes  of  units,  namely  cm, 
sec,  particle  (used  in  Reference  2)  and  meter,  sec,  kg  mole  (required 
for  the  present  analysis).  The  tables  are  self-explanatory,  though  a  word 
of  caution  is  in  order  for  the  ionization  reactions  10  through  15.  The 
rate  constants  for  these  reactions  sire  tabulated  as  the  forward  (i.e,  to 
the  right)  values,  which  is  contrary  to  the  remainder  of  the  appendix. 
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